SPECTROSCOPIC COSMOLOGICAL SURVEYS IN THE FAR-IR 
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Abstract 

We show the feasibility of spectroscopic cosmological 
surveys with the SAFARI instrument onboard of SPICA. 
The work is done through simulations that make use of 
both empirical methods, i.e. the use of observed luminosity 
functions and theoretical models for galaxy formation and 
evolution. The relations assumed between the line emis- 
sion to trace AGN and star formation activity have been 
derived from the observations of local samples of galaxies. 
The results converge to indicate the use of blind spec- 
troscopy with the SAFARI FTS at various resolutions to 
study galaxy evolution from the local to the distant (z~3) 
Universe. 

Specifically, two different and independent galaxy evo- 
lution models predict about 7-10 sources to be spectro- 
scopically detected in more than one line in a 2'x 2' SA- 
FARI field of view, down to the expected flux limits of 
SAFARI, with about 20% of sources to be detected at 
z>2. SPICA-SAFARI will be therefore excellent at de- 
tecting high-z sources and at assessing in a direct way 
their nature (e.g whether mainly AGN or Star Formation 
powered) thanks to blind spectroscopy. 

Key words: Galaxies: evolution, active galactic nuclei, star- 
burst - Missions: SPICA 



1. Introduction 

In the last years our perspective of galaxy evolution has 
greatly changed thanks to two main findings: A) the strong 
correlation observed in the local Universe between the 
mass of the black hole at the centre of a galaxy and the 
velocity dispersion of the stellar component of the galactic 
bulge (the so-called Magorrian relation)^ and B) the evi- 
dence that most, if not all, galaxies during their evolution 
pass through a FIR/submilhmetre bright phase. 



correlation for active galaxies in the local Universe, as the 
two processes appear - at least in local Seyfert galaxies 



(jTommasin et al L 12009 ) to be almost independent of each 
other. The Magorrian relation therefore has to find its 
origins at earlier epochs. 

The study of galaxy evolution necessarily implies in- 
vestigation of the full cosmic history of energy generation 
by stars (star formation and stellar evolution) and black 
holes (accretion), as well as of the energy loss processes 
such as the feedback from AGNs, because all these are 
responsible for the build up of the baryonic mass in the 
Universe and must ultimately be consistent and set up 
the observed local relation between luminosity and mass 
in galaxies. To understand the Magorrian relation we need 
to make the cosmic connections between stars in a galaxy 
and its massive black hole. The glo bal accretion power. 



measured at X-rays (|Hasinger et al L I2QQ5) and the star 



formation power, measured by Ra and rest- frame UV ob- 

times higher at 
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were 
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servations (jShim et al 
z=l-1.5 than today. On a cosmic scale, the evolution of 
supermassive black holes (SMBHs) appears related to the 



2004D 



ecestine: the presence of co-evolution ( Marconi et al . 


2004; 


Merloni et al . 


20041: Shankar et al 


12009). As sue-s-ested 


by various authors ffleckman et al . 


2004 


Granato et al . 



A) The Magorrian relation (Magorrian et al I . ll99^ 



2000l ) implies that the processes of black hole growth 
through mass accretion - and bulge formation - through 
star formation - are intimately linked. While large ellipti- 
cal galaxies with old stellar bulges are known to follow the 
Magorrian relation, we are not able to explain this tight 



the growth of bulges through SF may be directly 
linked to the growth of black holes through accretion. 
Quasars have been advocated as a source of negative feed- 
back that would quench star formation, however no clear 
evidence for this "negative feedback" has yet been found, 
while instead star formation is often, although not al- 
ways, concomitant with AGN/QSO activity over rather 
long duty cycles (100 Myr - few 100 Myr). It has also been 
realized that the optically bright phase of quasars covers 
not more than one tenth of their host lifetime, as obser- 
vations generally constrain quasar lifetimes to the range 
~ 10^ - 10^ yr. Optical studies of local galaxies show that 
FerraiM^fe: ^ffiPfi ff larg:e g;alaxy sphero ids host massive relic 

which, in turn, sug- 



black^holes (iKichstone et al . 199i 



gests that a mass-accreting AGN phase is one through 
which all galaxies pass. 

B) To account for the total energy generated by stars 
and that by black hole accretion, one first has to determine 
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how much of the observed luminosity is partly or heavily 2. MIR/FIR IMAGING SPECTROSCOPY IS THE KEY 
extinguished (reddened or obscured). We know from both 

galactic and extragalactic studies that the processes of star ^^^^ and young stars and black hole accretion disks show 
formation and early stellar evolution is usually deeply en- ^^^^o^g differences in the shape of their ionizing contin- 
shrouded in dust. Active galactic nuclei (AGNs) are often ^""i' However the far-UV continuum, dominating the to- 
opticallv obscured, both local ly e.g. (e.g. Goulding fc Alexan(!gn ^"l°"^^tric output luminosity in both processes, is in 
(3)) and at high redshifts (iHernan-Caballero et. al I . l2nn9h .g«"^^^l '^o* directly observable, due to absorption by HI. 
Thus, galaxies pass a significant fraction of their lifetime ^he best tracers and indeed discriminators of accretion 
deeply obscured by dust. 



Furthermore, recent observational studies assesses that 
during their evolution most (if not all) galaxies pass through 
a FIR/submillimetre bright phase, where star formation 
is probably the dominant energy production mechanism, 
but where also supermassive black- hole accretion may con- 
tribute significantly. The luminous and ultra-luminous in- 
frared galaxies, discovered in the 80's by the all-sky survey 
of IRAS, are rare in the local Universe, but they become 
much more common at high-redshift, as shown by the de- 
tection of the Submillimeter Galaxies (SMG). Their inte- 
grated luminosity accounts for a signifi cant fraction of th e 
cosmic infrared background (CIRB) (|Dole et al . 2QQ6h . 
the total energy of which exceeds that of the optical back- 
ground and implies that more than half the star-formation 
activity in the universe is hidden by dust extinction. 



and star formation processes are therefore emission lines 
from the photo-ionized gas. In both stabursts and AGN 
a fairly constant fraction (10-20%) of the ionising contin- 
uum gets absorbed by gas and then re-radiated as line 
emission, making this observable a powerful diagnostic 
tool. Detecting the exact fraction of the ionizing radiation 
absorbed by the gas surrounding the powering sources is 
however not crucial if one uses emission line ratios, which 
are therefore the most efficient observable probes of en- 
ergy production mechanisms. To overcome heavy extinc- 
tion, emission lines in the near-, mid- and far-IR have to be 
used, especially to probe obscured regions, such as those 
dominated by star formation and AGN activity. MIR/FIR 
imaging spectroscopy is therefore essential to trace galaxy 
evolution throughout cosmic times in an unbiased way and 
so to minimize dust extinction. 
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Figure 1. The positions of the infrared fine structure lines 
in the density -ionization diagram, showing the diagnos- 
tic potential of these lines to trace different astrophysical 
conditions: from photodissociation regions, to Stellar/HII 
regions, to AGN environments and high excitation coronal 
line regions fS p ino alio & Malkan . , .199k ). 
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Figure 2. Mid-IR spectra of Seyfe rt galaxies in the local 



Universe UTommasin et al , 200^). co mvared to those of 



starburst galaxies iBrandl et al 1 . \200di ). 



By the launch of SPICA, deep cosmological surveys 
undertaken by ISO, AKARI and Spitzer along with those 
planned for Herschel and SCUBA-2 will have produced 
catalogues containing the ffuxes of many tens of thou- 



sands of faint MIR/FIR/submm sources. Photometric sur- 
veys, however, can only be used to determine source counts 
and a great amount of multiwavelength photometric data 
and a priori assumptions are required for obtaining red- 
shifts and luminosities. Even more importantly, photomet- 
ric data alone do not allow to unambiguously differentiate 
between AGN and star formation activity. 

Substantial progress in studying galaxy evolution can 
only be achieved by using spectroscopic surveys. Having 
direct MIR/FIR spectroscopy can not only provide us with 
measured (rather than estimated) redshifts, but can also 
unambiguoulsly characterise the detected sources, measur- 
ing the AGN and starburst contributions to their bolomet- 
ric luminosities over a wide range of cosmological epochs. 
And this will all happen in a single observation. The great 
difference between the surveys that Herschel will perform 
and those made possible by an instrument like SAFARI 
is precisely this: the capability of obtaining unbiased spec- 
troscopic data in a rest-frame spectral range that has been 
shown in the local Universe: 1) to be unaffected by extinc- 
tion and 2) to contain strong unambiguous signatures of 
both AGN emission and star formation. 

Figure [1] shows how the IR fine structure lines cover 
the two main physical quantities - density and ioniza- 
tion - characterising photoionized and photon dissociated 
gas. Both star formation and stellar evolution processes 
as well as black-hole accretion can be traced from a de- 
voted study of a combination of lines and line ratios. The 
long wavelength of these lines ensures that dust extinction 
is minimized. S p ino glio fc M alka n (1992) first predicted 
line intensities of AGN and starburst galaxies using pho- 
toioization codes. To have an idea of the galaxies that 
could be observed by SAFARI in the far-IR in the red- 
shift range 0.4<z<4.0, one has to look at the rich rest 
frame mid-IR spectra of local galaxies that have been re- 
cently observed by the Spitzer satellite. Figure [2] shows 
the average Spitzer IRS high-resolution mid-IR spectra 
('Tommasin et al . 2QQ9I) of the 12/im sam ple of Seyfert 



galaxies (Rush. M alkan fc Spinoglio . 1992|). c ompared to 



992 

those of starburst galaxies ( Brandl et al L [2QQ6I ) . The qual- 
ity of the data is very high and shows the many features 
that can differentiate between AGN and star formation 
processes. 

3. Are local templates observable at high z ? 

The main question that we need to answer, before mov- 
ing to more detailed modeling, is if the local galaxies, for 
which we have detailed knowledge of their IR spectra, can 
be observed at high redshifts. In order to do so, we pre- 
dicted the line intensities as a function of redshift (in the 
range z=Q.l-5) fo r three local template ob jects: NGC1068 
(Alexander et al 1 l2nnnl : ISninoglio et all , ^b), a proto- 
typical Seyfert 2 gal axy containing bo th an AGN and a 
starburst; NGC6240 (|Lutz et al l . l2QQ3h. a bright starburst 
with an obscured AGN and M82 (|Forster Schreiber et al ' 



20011 : iGolbert et al I . ll999h . the prototypical starburst galaxy. 
We assumed that the line luminosities scale as the bolo- 
metric luminosity and adopted a luminosity evolution pro- 
portional to (z+1)^, consist ent with recent Spitzer lesui ts 



at least up to redshift z=2 (|Perez- Gonzalez et al I . l2QQ5f ). 
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Figures. Line observability with PACS onboard of Her- 
schel and with SPICA. Few selected diagnostic lines are 
shown as a function of redshift for the three template ob- 
jects M82, NGC 1068 and NGC6240 (from top to bottom). 
The lines have the same symbols as in Figure^ except for 
the open circles, which represent molecular lines. Line in- 
tensities are given in W m~'^ . The short dashed and long- 
dashed lines respectively give the 5a, 1 hour sensitivities 
of the PACS and SPICA spectrometers. 

For simplicit}0, we adopted an Einstein-De Sitter model 
Universe, with = and VtM= 1, IIo=75 km s~^ Mpc~-^. 
The luminosity distances have been derived using: 

dL{z) = {2c/Ho)-[l^z-(l^zY^^] 

. The results for the three template objects are reported in 
FigureO Among the brightest lines we show the [SIV]10.5/im, 
the [NeII]12.8/im and the [OIII]52/im diagnostic for the 
stellar/HII regions, the [NeV]24.3/im and the [OIV]25.9/im 
for the AGN component, the [OI]63/im and the [SiII]33.5/im, 



^ We note that the dependence on different cosmological 
models is not very strong. The popular model with Qm= 0.27, 
^A=0.73, Ho=71 km s~^ Mpc~^ would imply greater dilutions 
on line intensities, dilutions which increase with z by a factor 
of 1.5 for z=0.5 to 2.5 for z=5. 



for the photodissociation regions and the OH and H2 ro- 
tational hnes for the warm molecular component. 

It is clear from the figure that the PACS spectrome- 
ter will be able to observe only the most favorable object 
(NGC6240) up to z=2 in the brightest line ([OI]63/im), 
while the SPICA spectrometer goal sensitivities will al- 
low deep infrared spectroscopic studies of these relatively 
low luminosity objects (Lir=4, 20, 50 xlO^^L© for M82, 
NGC1068 and NGC6240, respectively) up to z - 1-2 for 
most lines and pushing this limit to even higher redshifts 
for the brightest lines. 

4. Starburst and AGN at z=0-4 with SAFARI 

With the goal of estimating the number of galaxies at 
different redshifts that can be detected by SAFARI spec- 
troscopic surveys, we have followed the method of: (1) 
identifying star formation and AGN activity by rest frame 
mid-IR spectroscopic tracers; (2) determining the correla- 
tions between line and continuum IR luminosity in the 
Local Universe; (3) using reliable IR continuum luminos- 
ity functions (derived from both observations and models) 
and (4) transforming the continuum luminosity functions 
(z=0-4) into line luminosity functions, from which to de- 
rive the number of detectable sources per redshift range. 
The use of different models (listed below) allows to provide 
a reliable range for SPICA- SAFARI predictions. 



Type 1 (Tommasin et al., 2009) 
Type 2 (Tommasin et al., 2009) 
Ulirg (Armus et al., 2007) 
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Figure 4. The [OIV]26iim/[SIII]33.5iim vs 

[Ne VI 24 . 3u m/fSIIIISS. 5u m diagram for local Seyfert 
galaxies l[ Tommasin et all \200A ) and ULIRG s 
( Armus etal\ . \20Q2)- 



4.1. Line ratio diagrams for AGN and starbursts 



A suitable line ratio diagram able to separate the emission 
from the AGN from that of the starburst can be built by 
using lines which can be produced only in the high ioniza- 
tion environment due to black-hole accretion (such as the 
[NeV] lines and, at a less extent, [OIV]26/im), together 
with lines generated by processes due to stellar forma- 
tion and evolution (such as those which are bright in HII 
regions: [Nell], [SIII], etc.). However, because we aim at 
covering a wide redshift range as possible with the SA- 
FARI spectrometer, we must choose lines which reside at 
long enough wavelengths. So we have concentrated on the 
[OIV]26/im/[SIII]33.5/im vs [NeV] 24.3/im/ [SIII] 33.5/im line 
diagnostic diagram, that we show in Figure [H because 
all these can be detected in the SAFARI spectral range 
for any redshift between z=0.4 and 5. In this figure we 
show that such line ratios can well separate AGN from 
star formation dominated galaxies, where we have plot- 
ted the line ratios of th e 12/im selected Seyfer^_galaxies_ 
( Tonimasin et al . 20091 ) and the ULIRG data ( Armus et al 
20071) . 



4.2. Correlations between line and continuum 

To be able to convert the IR and far-IR continuum lu- 
minosity functions into line luminosity functions, we have 
derived the correlations between line and continuum lu- 
minosities using the relatively large and complete sample 
of local Seyfert galaxies selected at 12/im which counts 
for more than 100 objects. We present in Figure [5] and 
Figure [6] the relation obtainedfor this sample between the 
far-IR luminosity vs the [SiII]34.8/im and [OIV]26/im line 
luminosities, respectively. From such a relation we then 
derived an analytical formula (obtained by least squares 
fits) to be used in the following analysis. 

4.3. Simulations with the Gruppioni model 
Our first method is based on the galaxy evolution model 



developed by Gruppioni and c ollaborators (jGruppioni fc Pozzi 



2OO9I : Gruppioni et al . 2009h . which makes use of all the 
available IR data to extrapolate continuum luminosity 
functions from z = to 2; ~ 4, with distinct contribu- 
tions from starbursts and AGN. To account for the popu- 
lation of spheroids (proto-ellipticals) which locally do not 
appear as substantial far-IR emitting sources, but are ex- 
pected to be relevant at high redshifts (e.g;. SMG g:a laxies), 
we integrated the model by Granato et al ('2004') which 
uses an ab-in itio approach for g alaxy formation and evolu- 
tion (see also iLapi et al I (|2006( )) capable to predict multi- 
wavelength luminosity functions and number counts. As 
an example, we show in Figure [71 the continuum and line 
luminosity functions and the predicted number of sources 
detectable in the [SiII]34.8/im line in the redshift range 




Figured. Correlation between the [SiII]34-8iim and 
the far-IR luminosity for the complete sample of 
12iim selected Seyfert galax ies in the local Universe 
iTommasin et alV \200^ . \200^ ). 



^12yLim Seyfert's (RMS93) 
log[L([0IV]26Aim)] = 1.062669658 * log [L(FIR)] 
R = 0.7659083594 
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Figured. Correlation between the [OIV]26fim and the far- 
IR luminosity for the complete sample of 12iim selected 
Sey fert galaxies in the local Universe jTommasin et al\ , 

\2m,Mm)- 



1.25<z<1.75, as derived from the modeled FIR luminos- 
ity function. 



Figure 7. Bottom pan el: Composite plot wi t h the FIR lu- 
minosity function from ICruppioni & Pozzi \200^ ) ( right 



hand side) and the derived [SillJS^-Sjim line luminosity 
function (left hand side) for the various galaxy populations 
in the redshift range 1.25<z<1.75: star forming galaxies, 
ULIG, Seyfert 2^s, Quasars and Spheroids. Top panel: 
derived total number of sources in a SAFARI field (2^ x 
2^) at 1.25<z<1.75 as a function of the sensitivity. The 
vertical line shows the 5a, 1 hour sensitivity limit of 2.0 
X 10-^^ W m-^. 



By following this approach, we expect to detect about 
7 sources per field of view, over the redshift range 0.<z<4.0, 
in the [SiII]34.8/im line, which is a star formation indi- 
cator. We assumed a 5cr, 1 hour detection limit of 2x 
10-19 w The lines of [OIV]26/im (indicator of AGN 
activity) and [SIII]33/im are, in the local galaxy popula- 
tion, a factor 1.5-2 fainter than the [Sill] line, though they 
are found to b e of the same order o f magnitude in LIRGs 
and ULIRGs ( Veilleux et al . 2QQ9I ). Given the above, we 
expect to detect all three lines simultaneously in all the 
observed sources. 

4.4. Simulations with the Franceschini model 

As a second approach we adopted the model developed 
by Franc eschini an d collaborators at Padova University 
(iFranceschini et al 112009). This is a backward evolution 
model which fits all available data from Spitzer, ISO, COBE, 
SCUBA, etc. It includes direct determinations of multi- 
wavelength redshift-dependent luminosity functions from 
Spitzer and accounts in great detail not only for star form- 
ing galaxies, but also for type-1 and type-2 AGNs. 
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Figure 8. Prediction of the results of a spectroscopic sur- 
vey over 0.5 square degrees (requiring of order of 500 hrs 
integration with n o mina l sensitivity) from the model by 
Franceschini et al \200^ ). The number of detected sources 
would be 120 Type-1 AGNs, 770 Type-2 AGNs and 1870 
starburst galaxies in Sill, for a total of 2800 objects. It 
has been adopted the goal sensitivity of SAFARI of 2.0 x 
W m-'^ (5a, 1 hour). 

The simulation has considered the contributions of: 
moderate-luminosity starburst galaxies and high-luminosity 
starbursts, for which we have considered the [SiII]34.81/im 
line; Type-1 Active Galactic Nuclei and quasars, for which 
we have considered the [NeV]24.3/im and [OIV]26/im lines. 
Type-2 AGNs also provide substantial contributions to the 
detection statistics (they are considered to be a fraction of 
the low-L and high-L SBs) We referred to the 2.0 x 10~^^ 
W m~^ 5<j, Ihr limit. We predict ~4 starburst, 2 type-2 
and 1 type-1 AGN per SAFARI field of view of 2'x 2( for 
a total of 7 sources, in good agreement with the predic- 
tions obtained in Section 4.3 with a totally independent 
model. We show, as an example, in Figure [8] the results of 
a spectroscopic survey over 0.5 square degrees. 

5. Conclusions 

The main conclusions of this work are: 

- After decades of efforts we are close to having reli- 
able IR measures of star formation rate and AGN ac- 
cretion power ^ through IR/FIR spectroscopic surveys^ 
completely unaffected by dust, allowing us to study the 
evolution of galaxies in terms of the power produced 
by the main energy production mechanisms. 

- Accurately measuring the power due to stellar forma- 
tion & evolution and that due to gravity in collapsed 
nuclei is the first step towards understanding them, 
and how they have been related over the history of the 
Universe. 



- We learned how to measure these in local galaxies 
through mid-IR spectroscopy. 

- Blind FIR spectroscopic surveys with SAFARI can be 
the mean to physically measure galaxy evolution, pro- 
vided that a sensitivity of 2x 10~^^ W/m^ (5cr, 1 hour) 
can be reached. 

- Two different and independent galaxy evolution mod- 
els predict about 7-10 sources to be spectroscopically 
detected in more than one line in a 2' x 2' SAFARI field 
of view, down to the expected fiux limits of SPICA, 
with about 20% of sources to be detected at z>2. 
SPICA-SAFARI will be therefore excellent at detect- 
ing high-z sources and at assessing in a direct way their 
nature (e.g whether mainly AGN or Star Formation 
powered) thanks to blind spectroscopy. 
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